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SUMMARY  PACE 


THE  PROBLEM 

The  variety  of  stimulus  waveforms  used  in  vestibular  research  and 
in  clinical  testing  requires  a  versatile  procedure  for  quantitatively 
analysing  oculomotor  responses. 

FINDINGS 

A  new  procedure  was  developed  to  analyse  oculomotor  responses  to 
many  forms  of  vestibular  and/or  visual  stimuli  used  in  experimental  work 
and  in  clinical  practice.  The  main  advantage  of  the  procedure  lies  in 
the  simplicity  of  the  basic  algorithm  which  enables  fairly  rapid  data 
analysis  of  a  variety  of  response  forms  on  a  medium-speed  computer 
working  in  a  high-ievel  language.  For  slow  phase  analysis,  detection  of 
saccades  is  carried  out  by  a  simple  threshold  procedure  based  upon  the 
expected  response  form,  and  regression  procedures  are  used  to  obtain 
stimulus-response  correlation  and  slope  (gain)  measures.  The  procedure 
provides  measures  of  gain,  phase  and  directional  preponderance  for 
responses  to  sinusoidal  stimuli. 
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This  report  was  written  while  the  author  was  a  visiting  scientist  at  the 
Naval  Aerospace  Medica.i  Research  Laboratory.  The  algorithm  was  adapted 
by  Dr.  Barnes  to  NAMRL  equipment  which  is  similar,  though  not  identical, 
to  data  processing  equipment  at  the  Royal  Mr  Force  Institute  of  Aviation 
Medicine,  Farnborough,  England.  The  major  work  on  program  development 
was  accomplished  by  Dr.  Barnes  before  his  visit. 
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INTRODUCTION 


There  are  already  many  examples  of  procedures  for  the  analysis  of 
nystagmus  (1,2,13),  some  of  which  are  able  to  perform  on-line  analysis 
with  a  suitably  fast  computer.  Most  of  these  procedures  involve  the  use 
of  some  fairly  complex  algorithm  for  distinguishing  becween  fast  and 
slow  phases  of  the  eye  movement.  The  motivation  for  developing  a  differ¬ 
ent  procedure  in  our  laboratory  was  threefold.  The  first,  and  foremost, 
reason  was  an  essentially  practical  one.  As  a  research  group  with 
continually  varying  computational  needs,  we  opted  for  the  installation 
of  a  fairly  powerful  desktop  computer  (Hewlett-Packard  9845S)  which 
could  be  easily  programmed  by  experimenters  in  a  high-level  language 
(BASIC)  without  the  assistance  of  a  dedicated  computer  programmer.  The 
drawback  inherent  in  the  use  of  such  a  computer  is  that  it  is  relatively 
slow  compared  with  those  which  use  machine-code  instructions.  There  was 
thus  an  incentive  to  develop  an  analysis  procedure  which  would  combine 
the  superior  pattern-recognition  characteristics  of  the  human  operator 
with  a  relatively  simple  method  of  detecting  the  occurrence  of  fast 
phases . 

The  second  reason  for  developing  a  new  procedure  came  about  as  a 
result  of  the  need  to  analyse  the  oculomotor  responses  to  voluntary  head 
movements.  In  these  responses  the  peak  slow-phase  eye  velocity  was 
frequently  very  high  (>  100  deg/sec)  and  the  normal  procedures  for 
detecting  fast  phases,  which  might  frequently  be  reduced  to  100-150 
deg/sec,  were  inadequate.  The  third  point  was  that  a  reliable  method 
was  required  for  detecting  fast  phases  in  the  oculomotor  response  to 
many  different  types  of  stimuli.  In  particular,  there  was  a  need  to  be 
able  to  analyse  the  responses  to  a  high  frequency  visual  pursuit  task  in 
which  the  fast  phases  are  likely  to  be  in  the  same  direction  as  the 
slov.’-phase  eye  movements. 

There  are  two  means  by  which  the  improvement  in  processing  time  has 
been  achieved.  The  first  is  that  the  detection  of  saccades  is  carried 
out  by  a  simnle  threshold  procedure  based  upon  the  expected  form  of  the 
response.  The  threshold  limits  can  be  set  up  by  the  operator,  although 
an  automatic  procedure  is  available  for  those  responses  which  are  unlikely 
to  differ  greatly  from  the  population  norm.  The  second  feature  is  that 
there  is  no  attempt  in  the  analysis  to  interpolate  between  slow  phases 
in  order  to  reconstitute  the  cumulative  slow-phase  eye  position.  Instead, 
the  points  corresponding  to  the  fast-phase  components  are  simply  discarded, 
leaving  only  slow-phase  eye  velocity  points  for  subsequent  statistical 
analysis . 

THE  FAST- PHASE  DETECTION  PROCEDURE 

DATA  PREPARATION 

Prior  to  any  data  analysis  a  number  of  analogue  and  digital  techniques 
arc  used  to  prepare  the  data.  First:  the  eye  movement  recording  signal 
(normally  an  el  o.c  t  ro-ocu  1  ograph  i c  recording)  is  filtered  by  an  analogue 
.low-pass  filter  (Kento  VBF/19)  with  a  48  dli/oct  ave  roll-off  at  50  Hz,  If 
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the  stimulus  is  also  to  be  recorded,  as  it  is  during  an  oscillatory 
stimulus,  then  it  is  also  filtered  in  a  similar  manner  so  that  any  phase 
distortion  of  the  waveform  will  apply  equally  to  both  stimulus  and 
response.  The  filtered  signals  are  passed  through  an  analogue-digital 
converter  (Microconsultants  Modular  15)  which  is  capable  of  sampling  up 
to  eight  channels  simultaneously,  using  a  parallel  sample  and  hold 
module.  The  sampling  rate  and  number  of  channels  to  be  sampled  can  be 
set  up  prior  to  data  acquisition  under  computer  control. 

The  digitised  data  are  feu  directly  into  the  computer  (Hewlett- 
Packard  9845S)  using  Direct  Memory  Access,  a  rapid  procedure  with  a 
transfer  rate  of  20,000  words/sec.*  Analogue  signals  obtained  from 
experiments  are  normally  recorded  on  a  multichannel  tape  recorder  rather 
than  being  digitised  on-line  and  stored  in  the  computer.  It  is  felt 
that  this  method  gives  greater  flexibility  should  it  ever  be  necessary 
to  analyse  the  records  in  a  different  manner.  The  advantage  of  using 
Direct  Memory  Access  is  that  it  enables  the  analogue  tape  to  be  replayed 
at  many  times  the  recording  rate  (up  to  32  x  on  our  own  recorder  -  a 
Racal-Thermionic  Store  4-D)  and  so  drastically  reduce  the  data  '’cquisition 
t  itne . 

If  the  recorded  eye  movement  signal  is  sti’l  noisy  despite  the 
analogue  filtering,  the  computer  program  incorporates  a  subroutine  to 
apply  a  digital  filter  as  an  operator-specified  option.  The  filter  is 
of  the  auto-regressive  type  developed  by  Lynn  (10,11).  The  filter 
introduces  no  phase  distortion,  but  great  care  must  be  used  in  selecting 
the  original  sampling  rate  if  this  procedure  is  "o  be  used.  It  is 
particularly  useful  when  attempting  to  analyse  EOG  signals  Induced  by  a 
periodic  stimulus  of  frequency  greater  than  1  to  2  Hz,  where  the  signal 
to  noise  ratio  may  be  severely  diminished. 

All  the  operations  performed  on  the  eye  movement  signals  are  actually 
carried  out  on  the  eye  velocity  signal.  Thus,  before  any  further  proces¬ 
sing  can  take  place  the  eye  displacement  signal  is  digitally  differentiated, 
using  a  two-point  difference  method,  as  shown  below: 

0  i  =  (0i  +  i  “  0i  _  i}  7  2-Dt 


*Equipment  components  at  NAMRL  which  differ  from  some  of  those  described 
in  this  section  were  adapted  by  the  author  to  accomplish  versatile  quanti¬ 
tative  eye  movement  analysis.  At  NAMRL  the  analog  signals  are  recorded  on 
a  cassette  recorder  (Phillips  Mini-Log  4  or  a  TEAC  R  80).  These  signals, 
without  analog  filtering,  are  played  back  into  a  Hewlett-Packard  multi- 
programmer,  Model  6942A,  which  incorporates  three  plug-in  cards:  an  analog- 
to-digital  converter,  timer-pacer,  and  memory  card.  The  multiprogrammer 
under  the  control  of  the  Hewlett-Packard  9845T  computer  samples,  digitizes, 
and  stores  the  data  on  the  memory  card  where  they  are  transferred  to  the 
computer  memory  upon  request.  This  process  is  limited  by  memory  card 
capacity  to  4096  readings  per  channel,  and  is  slower  than  the  direct  memory 
access  method  used  at  IAM. 
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where  0^  =  eye  velocity  at  time  t  = 

0^  =  eye  displacement  at  time  t  =  i.Dt 
Dt  -  sampling  increment(s) 

If  the  stimulus  is  a  periodic  waveform,  then  it  is  also  differentiated 
and  stored  along  with  the  original  waveform  to  be  used  in  the  subsequent 
statistical  analysis. 

The  eye  velocity  signal  is  presented  on  the  computer  graphics 
display  as  shown  in  Figure  1.  In  this  example,  in  which  the  stimulus  is 
a  sinusoidal  waveform  at  0.1  Hz,  the  stimulus  velocity  is  also  displayed. 
This  exhibits  the  typical  response  of  the  vestibulo-ocular  reflex,  in 
which  the  slow-phase  eye  velocity  is  modulated  in  accord  with  the  stimu¬ 
lus  waveform  but  is  in  the  compensatory  direction.  (N.li.  The  polarity 
of  the  stimulus  waveform  is  reversed  in  Figure  1.)  The  slow  phase  is 
interspersed  with  brief  high  velocity  spikes  of  reversed  polarity, 
representing  the  differentiated  fast  phases.  Figure  1  also  shows,  for 
completeness,  the  eye  displacement  signal,  though  it  is  not  normally 
displayed  to  the  operator  during  the  analysis  on  the  computer  visual 
display  unit  (VDU). 

SETTING-UP  THRESHOLD  LIMITS  FOR  SINUSOIDAL  STIMULI 

The  procedure  which  is  used  for  fast-pha3e  detection  is  one  in 
which  the  expected  response  is  used  to  set  up  velocity  thresholds  for 
the  assessment  of  slow  and  fast  phases.  The  procedure  is  illustrated  in 
Figure  1  for  a  sinusoidal  stimulus,  although  the  procedure  is  similar 
for  other  types  of  stimuli  as  will  be  shown  later. 

The  first  step  is  to  estimate  the  characteristics  of  the  stimulus 
waveform.  The  frequency  is  calculated  by  finding  the  zero  cross-over 
points,  and  the  peak  velocity  is  estimated  by  integration  of  the  first 
half-cycle  of  the  stimulus  waveform.  Only  a  complete  number  of  cycles 
is  used  for  the  analysis,  which  enables  any  offset  in  the  stimulus 
waveform  signal  to  be  identified  and  removed.  The  graphics  display  has 
a  cursor  facility  in  which  the  position  of  one  horizontal  and  one  vertical 
line  can  be  made  to  move  in  a  direction  perpendicular  to  the  line.  In 
the  manual  mode  of  threshold  determination  these  cursor  lines  are  brought 
up  on  the  display  screen  so  as  to  overlay  the  eye  velocity  waveform. 

The  first  set  is  aligned  with  the  vertical  cursor  indicating  the  position, 
of  the  first  peak  of  the  stimulus  waveform  and  with  the  horizontal 
cursor  indicating  the  estimated  level  of  peak  slow-phase  eye  velocity. 

This  estimate  is  taken  to  be  the  peak  stimulus  velocity  level  multiplied 
by  the  normal  gain  factor  (approximately  0.7)  for  the  vestibulo-ocular 
reflex  when  recorded  in  the  dark  (fi.  If  the  position  of  the  vertical 
cursor  does  not  coincide  with  the  peak  slow-phase  eye  velocity,  or  in 
other  words,  if  there  is  a  phase  difference  between  stimulus  and  response, 
the  cursor  is  adjusted  appropriately.  Likewise  if  the  horizontal 
cursor  does  not  give  a  good  estimate  of  peak  eye  velocity,  it  can  be 
maneuvered  in  order  to  do  so.  A  second  set  of  cursors  then  allows  the 
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Figure  1 

An  illustration  of  the  manner  in  which  stimulus-related  threshold 
limits  are  set  up  to  dissociate  the  fast-  and  s low-phase  components 
of  nystagmus  elicited  in  response  to  a  0.1  Hz  sinusoidal  oscillation 
of  the  head  about  the  yaw  axis. 

(a)  Eye  position  signal  (EOG)  sampled  at  20-ms  intervals.  (N.B.  The 
steplike  nature  of  the  signal  is  a  common  aberration  produced  by 
such  computer  graphics  systems.) 

■b)  The  differentiated  EOG  signal  with  fast  and  slow  phases  present. 
U)  The  eye  velocity  signal  after  fast  phases  have  been  removed. 

(d)  The  stimulus,  reversed  in  polarity  for  simplification  of  analysis. 
The  analysis  is  normally  carried  out  on  ten  cycles  of  the  response, 
although  only  three  are  displayed  here  for  clarity. 
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position  and  magnitude  of  the  troughs  of  the  response  to  be  estimated. 

If  there  is  my  directional  preponderance  in  the  response,  this  method 
allows  the  difference  between  the  peak  and  trough  velocities  to  be  used 
as  an  estimate  of  that  offset. 

In  this  manner  an  estimate  of  peak  slow-phase  eye  velocity  together 
with  any  offset  and  phase  change  with  respect  to  the  stimulus  waveform 
have  been  obtained.  The  next  step  is  to  use  this  estimate  to  set  up 
threshold  limits  on  either  side  of  the  estimated  response  waveform.  This 
threshold  is  normally  taken  to  be  +  50  deg/sec,  and  a  set  of  small 
cursor  lines  is  brought  up  to  indicate  where  the  threshold  line  will 
pass  in  relation  to  the  first  peak.  If  this  threshold  level  is  felt  to 
be  insufficient  for  the  noise  level  of  the  eye  movement  signal,  then  it 
can  be  increased  by  moving  the  cursor  the  appropriate  distance  in  the 
direction  of  the  ordinate.  The  threshold  is  not  normally  decreased 
below  50  deg/sec  when  analysing  eye  movements  recorded  by  EOG  because  of 
the  inherent  noise  in  the  signal.  When  a  suitable  threshold  hrs  been 
assigned,  the  threshold  limits  are  plotted  by  the  computer  graphics  in 
relation  to  the  eye  velocity  waveform,  as  shown  in  Figure  1. 

THRESHOLD  LIMITS  FOR  TRANSIENT  RESPONSES 

In  many  types  of  vestibular  stimuli  the  oculomotor  response  is 
transient,  decaying,  or  changing  with  respect  to  time  in  a  manner  which 
is  not  always  easy  to  correlate  with  the  stimulus.  In  such  examples, 
two  methods  of  setting  threshold  limits  have  been  adopted.  First,  for 
those  responses  which  are  reasonably  predictable,  an  estimate  of  the 
slow-phase  velocity  is  obtained,  based  upon  established  population 
norms.  For  example,  the  response  of  the  lateral  semicircular  canals  to 
a  rapid  deceleration  from  a  constant  angular  velocity  exhibits  an 
exponent  i.al  decay  of  slow-phase  eye  velocity  which  is  governed  pre¬ 
dominantly  by  a  time  constant  of  approximately  15  seconds,  as  shown  in 
Figure  2.  Accordingly,  the  threshold  limits  may  be  set  up  on  the  assump¬ 
tion  that  tlie  response  exhibits  such  a  form  (Figure  2),  and  this  will  be 
adequate  even  for  quite  large  variations  in  the  actual  time  constant. 
Second,  if  tlie  response  is  not  easily  predicted  in  terms  of  some  simple 
mathematical  representation,  the  operator  is  allowed  to  manipulate 
cursors  in  order  to  simulate  tlie  slow-phase  velocity  with  as  many  strnight- 
.1  Ine  segments  as  are.  necessary.  This  technique  is  particularly  useful 
in  analysing  the  responses  to  a  caloric  irrigation  as  shown  in  Figure  3. 

THE  ELIMINATION  OF  FAST  PHASES 

Every  point  which  lies  outside  the  threshold  limits  is  taken  to 
represent  a  fast-phase  eye  movement.  The  advantage  of  this  procedure  is 
that  Lt  allows  fast  phases  of  fairly  low  velocity  (s  50  deg/sec)  to  be 
detected.  One  of  the  problems  inherent  in  my  procedure  which  operates 
on  digitised  data  is  that  the  process  of  sampling  the  data  effectively 
filters  it.  Thus  the  peak  level  of  the  digitised  data  during  a  fast 
phase  may  not  represent  true  peak.  This  is  illustrated  in  Figure  4, 
where  sampling  point  8  is  seen  not  to  occur  at  actual  peak  velocity. 

This  effect  is  made  worse  when  the  sampling  rate  is  reduced.  Most  last 
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phases  in  vestibular  nystagmus  last  for  20-60  ms  (14).  It  is  therefore 
necessary,  in  normal  circumstances,  that  the  sampling  interval  not  be 
increased  beyond  20  ms  and  it  should  preferably  be  10  ms  or  less.  If 
any  analysis  of  fast-phase  characteristics  is  to  be  carried  out,  the 
sampling  interval  should  be  approximately  5  ms  or  less.  However,  even 
if  there  are  only  one  or  two  samples  per  fast-phase  movement,  the  tech¬ 
nique  for  threshold  detection  still  works  adequately. 

The  sampled  points  (e.g.,  7,  8,  9,  and  10  in  Figure  4)  that  lie 
outside  the  threshold  limits  are  automatically  discarded  from  the  analy¬ 
sis,  In  addition,  in  order  to  ensure  that  no  component  of  the  fast 
phase  remains,  one  point  preceding  the  first  point  above  threshold 
(point  6  in  Figure  4).  and  one  point  following  the  last  point  above 
threshold  (point  11  in  Figure  4)  are  also  discarded.  The  rate  at  which 
the  signal  is  sampled  is  critical  in  this  process.  It  has  been  deter¬ 
mined  empirically  that  the  number  of  occasions  on  which  there  is  any 
remnant  of  the  fast  phase  is  negligible  if  the  sampling  interval  is 
greater  than  or  equal  to  20  ms  if  the  threshold  is  set  at  50  deg/sec. 

This  is  partially  supported  by  the  recorded  velocity  profiles  of  sac¬ 
cadic  eye  movements  (8)  in  which  the  velocity  tends  to  exceed  50  deg/sec 
well  within  the  first  10  ms  of  the  response.  Consequently,  if  the 
sampling  interval  is  less  than  10  ms,  the  number  of  points  below  the 
threshold  which  are  discarded  should  be  increased  appropriately.  Con¬ 
versely,  if  the  sampling  interval  is  too  great  (>  40  ms),  the  process  of 
casting  out  data  points  may  result  in  some  slow-phase  points  being 
unnecessarily  lost. 

This  procedure  for  eliminating  the  fast  phases  is  common  to  all  the 
methods  of  analysis.  However,  the  reduced  data  are  statistically  analysed 
in  different  ways  for  different  types  of  stimuli.  In  the  following 
section,  an  outline  of  these  procedures  is  given. 

ANALYSIS  OF  REDUCED  DATA 

ANALYSIS  OF  SINUSOIDAL  RESPONSES  USING  REGRESSION  PROCEDURES 

One  of  the  problems  associated  with  discarding  data  in  the  manner 
just  described  is  that  the  data  cannot  be  subjected  to  the  more  conven¬ 
tional  forms  of  frequency  domain  analysis,  such  as  the  Fourier  Transform, 
because  the  data  points  are  no  longer  equi-spaced.  However,  it  is  still 
possible  to  obtain  a  direct  correlation  between  the  stimulus  and  the 
response  using  a  regression  procedure.  (This  may  be  best  illustrated,  as 
in  Figure  5,  by  plotting  eye  velocity  (0)  against  head  velocity  (I)  for 
data  of  the  type  exhibited  in  Figure  1.  If  the  eye  velocity  is  almost 
exactly  compensatory  for  head  velocity  (l.e,,  there  is  negligible  phase 
error),  then  the  slow-phase  eye  movements  will  be  disposed  in  a  straight 
line  which  passes  through  the  origin.  The  slope  of  this  line,  which 
should  be  negative,  represents  the  gain  of  t Lie  response.  In  Figure  5, 
the  threshold  limits,  defined  in  Figure  1,  appear  as  straight  lines 
enclosing  the  slow-phase  data  points.  The  fast-phase  components  lie 
outside  the  threshold  within  the  circled  areas. 
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Figure  4 

An  illustration  of  the  method  by  which  fast-phase  eye  movements  are 
detected  and  discarded,  using  simple  threshold  limits  on  the  eye  velocity 
waveform.  When  a  data  po'nt  (e.g.,  point  7)  lies  outside  the  threshold, 
both  that  point  and  points  within  the  preceding  20  ms  are  discarded  from 
the  analysis. 
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Since  eye  velocity  (0)  should  be 'proportional  to  head  velocity  ($) , 
it  is  possible  to  carry  out  a  linear  regression  on  the  slow- phase  points 
to  obtain  an  accur  te  determination  of  the  slope  and  intercept.  For 
this  purpose  the  sampled  data  points  of  the  stimulus  waveform  which  are 
paired  with  the  fast  phase  of  the  eye  velocity  waveform  ^re  also  discarded. 
However,  if  there  .is  any  phase  difference  between  0  and  <£ ,  the  slow- 
phase  points  in  Tigure  5  appear  as  an  ellipse,  not  a  straight  line.  It 
is  possible  co  obtain  an  indication  of  this  phase  error  if  the  correlation 
is  made  between  eye  velocity,  head  velocity,  and  head  acceleration  ^ . 

That  is,  it  is  assumed  that  the  data  points  will  fit  an  equation  of  the 
form: 

0  -  A<j>  +  B<j>  +  C 

where  A,  B,  and  C  are  constants  of  proportionality. 

If  head  velocity  is  a  sinusoidal  function  then 

«  Ik 

<t>  =  sin  mt  and  <j>  =  ftw  cos  rnt  (2) 

where  SI  *  peak  stimulus  velocity  (deg/sec). 

Substituting  these  values  in  equation  (1)  gives 

0  =  ASi  sin  uit  +  BQrn  cos  ut  +  C  (3) 

=  \/a2  +  B2  u)2'  SI  sin  (cut  -  ip)  +  C  (4) 


where  t|j  =  tan 


-1  Bw 


i/2  2  2’ 

In  this  manner  it  is  possible  to  extract  the  gain(\A  +  B  w  ), phase  (ty), 
and  directional  preponderance  (C)  from  the  response.  In  addition,  it  is 
possible  to  apply  more  complex  regression  equations  to  extract  any 
nonlinearities  from  the  response. 

ANALYSIS  OF  TRANSIENT  RESPONSES  USING  MODEL  PREDICTIONS 

As  discussed  earlier,  some  stimuli  applied  to  the  vestibular  apparatus 
give  rise  to  transient  eye  movement  responses  which  cannot  be  directly 
correlated  with  the  stimulus.  However,  it  is  possible  to  establish 
models  of  oculomotor  performance  (5)  which  have  a  sound  physiological 
basis,  and  the  response  can  be  fitted,  using  regression  procedures,  to 
such  a  model.  The  models  used  may  be  much  more  complex  than  those 
simple  representations  which  were  first  used  to  establish  threshold 
limits . 


An  example  of  the  type  of  response  which  can  be  fitted  in  this 
manner  is  that  shown  in  Figure  6.  The  stimulus  is  a  ramp-step  angular 
velocity  function  (an  acceleration  pulse)  applied  to  the  horizontal 
semicircular  canals.  The  slow-phase  eye  velocity  rises  to  a  peak  value 
at  the  end  of  the  acceleration  and  then  decays  in  an  exponential  manner. 
The  response  cannot  be  modeled  by  a  single  exponential  since  the  reversal 
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after  20-30  seconds  indicates  the  presence  c  a  second  time-dependent 
component  which  is  normally  referred  to  as  adaptation  (6,15).  The 
response  can  be  adequately  modeled  (5)  using  a  transfer  function  of  the 
form 

9  -  KGTBTAD_S  0 

V  (1+Tbs)(1+Tads) 

where  K  is  a  constant  of  proportionality  (approx.  0.7) 

c 

T  is  the  time  constant  of  adaptation  (approx.  80-100 
AD  N 
sec) 

T  is  the  long  time  constant  of  the  canal-ocular 
^  response  (approx.  15  sec) . 

The  response  of  equation  (5)  to  an  acceleration  pulse  (level  a)  is 

6  =  a  Kc  F(t) 


where  F(t)  ■  (e 


(1  -  e 


'■  for  0  <  t  < 

AD  (  t/TB.  “t/TB 
-  (1  -  e  )  e 


t  <  t 


This  equation  can  be  fitted  to  the  data  by  assuming  initial  estimates 
for  l’AD  and  Tfi  and  then  using  an  iterative  multiple  regression  procedure 
(3)  to  obtain  new  estimates  of  T  ,  T  ,  and  K  which  are  within  acceptable 
limits  of  accuracy. 

This  procedure  c.-n  be  carried  out  without  modification  on  the  data 
points  which  remain  after  the  fast-phase  components  have  been  discarded. 
However,  this  frequently  involves  a  large  number  of  data  points  which 
can  result  in  a  long  processing  time  because  of  the  iterative  procedure. 
Consequently,  a  technique  for  reducing  the  data  points  is  normally  used. 
This  is  carried  out  by  averaging  the  number  of  data  values  which  are 
remaining  within  consecutive  time  bins  of  specific  duration  (say,  1  sec 
for  the  type  of  response  shown  in  Figure  6) ,  a  method  which  is  only 
suitable  for  responses  varying  slowly  with  time. 

ANALYSIS  OF  TRANSIENT  RESPONSES  USING  ARBITRARY  DESCRIBING  FUNCTIONS 

On  occasions  it  is  only  necessary  to  extract  one  or  two  important 
features  from  a  response  rather  than  to  characterise  it  in  terms  of  a 
meaningful  model.  An  example  is  the  caloric  response  in  which  the  peak 
slow-phase  velocity  and  the  time-to-reach  peak  are  two  of  the  most 
important  variables  which  can  be  used  to  identify  abnormalities  of 
vestibular  function.  Models  of  the  caloric  response  have  been  developed 
(6,12)  but  are  complex  and  would  not  be  easily  fitted  using  a  regression 
procedure.  However,  in  order  to  extract  the  important  variables  it  is 
sufficient  to  fit  the  data  with  a  simple  power  law  regression.  A  5th 
order  solution  is  normally  used  for  caloric  analysis,  giving  a  fit 
through  the  data  points  as  shown  in  Figure  7.  Peak  slow-phase  velocity 


Stow  phase 
eye  velocity 


The  s low-phase  eye  velocitj  obtained  by  a  caloric  irrigati 
power  function  of  time,  from  which  peak  velocity  and  the  t 


and  time  to  peak  are  then  easily  derived  from  the  parameters  of  the 
fitted  equation. 


CONCLUSIONS 

A  new  procedure  for  analysing  nystagmus  and  other  oculomotor  responses 
has  been  described.  The  main  advantage  of  this  procedure  lies  in  the 
simplicity  of  the  basic  algorithm  which  enables  the  data  to  be  analysed 
fairly  rapidly  on  a  medium-speed  computer  working  in  a  high-level 
language.  The  method  has  been  found  particularly  useful  for  analysing 
responses  of  the  vestibulo-ocular  reflex  in  which  the  slow-phase  eye 
velocity  may  be  very  high,  as,  for  example,  during  voluntary  head  move¬ 
ments,  but  is  also  applicable  to  many  other  forms  of  vestibular  responses, 
to  pursuit  tracking  responses,  and  to  optokinetic  responses.  In  addi¬ 
tion,  because  the  analysis  is  carried  out  on  eye  velocity  data  and 
because  regression  procedures  are  used  for  statistical  analysis,  it  is 
not  necessary  to  manufacture  new  data  points  to  interpolate  across  the 
gaps  left  by  the  discarded  fast  phases. 

The  principal  disadvantage  of  the  technique  described  here  is  that 
unless  a  reasonable  prior  estimate  of  the  slow-phase  eye  velocity  can  be 
made,  it  is  unsuitable  for  on-line  analysis.  This  problem  can  be  over¬ 
come  for  certain  types  of  stimuli  such  as  caloric  irrigation  in  which 
the  slow-phase  response  is  normally  only  in  one  direction  and  is  of 
relatively  low  velocity  compared  with  the  fast  phases.  It  has  also  been 
found  that  in  experiments  on  normal  human  subjects,  the  response  differs 
sufficiently  little  from  that  of  the  population  norm  that  analyses  can 
be  carried  out  automatically.  One  further  aspect  of  this  procedure 
which  carries  a  slight  penalty  is  that  in  order  to  cut  down  the  amount 
of  data  points  stored  and  thus  decrease  processing  time,  the  sampling 
rate  must  be  kept  as  low  as  possible  (optimally  10-20  ms).  This  precludes 
the  possibility  of  analysing  the  characteristics  of  the  fast-phase  eye 
movements  in  any  detail  and  also  makes  it  difficult  to  analyse  nystagmus 
in  which  the  saccadic  frequency  is  very  high  (e.g.,  in  older  patients 
an',  those  using  nicotine).  However,  these  two  problems  may  still  be 
adequately  addressed  if  a  longer  processing  time  is  tolerable. 

The  procedure  described  here  has  been  used  to  analyse  the  responses 
in  many  different  forms  of  experiments  (4,7,9)  and  in  routine  clinical 
analyses  and  has  been  found  to  be  reliable  in  its  detection  and  elimina¬ 
tion  of  fast  phases.  Comparisons  of  analyses  carried  out  in  this  manner 
with  those  carried  out  by  hand  using  a  beat-by-beat  slope-fitting  procedure 
have  established  that  the  method  is  accurate.  In  addition,  it  is  less 
susceptible  to  the  inclusion  (or  exclusion)  of  data  values  which  the 
human  operator  subjectively  might  feel  should  (or  should  not)  be  included 
in  the  analysis. 
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